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Abstract 

We examine in detail the possibilty of applying Darboux transformation to non Hermitian hamiltoni- 
ans. In particular we propose a simple method of constructing exactly solvable VT symmetric potentials 
by applying Darboux transformation to higher states of an exactly solvable VT symmetric potential. It is 
shown that the resulting hamiltonian and the original one are pseudo supersymmetric partners. We also 
discuss application of Darboux transformation to hamiltonians with spontaneously broken VT symmetry. 
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1 Introduction 



Ever since it was conjectured by Bender et al that some non-Hcrmitian Hamiltonians exhibiting symmetry 
under the combined transformation of parity (V : x — > —x), and time reversal (T : i —> —i) admit real 
eigenvalues non-Hermitian Hamiltonians have been the basis of many recent works on VT symmetry 

and pseudo-Hermiticity |3lE]j because of intrinsic interest and their possible applications in molecular 
physics, quantum chemistry, superconductivity, quantum field theory and others. 

On the other hand there are not many examples of exactly solvable complex potentials (both VT 
invariant as well as otherwise). However, as in the Hermitian case, there have been attempts to expand 
the class of exactly solvable non Hermitian potentials by using different methods [SUSIE]. In the Hermitian 
case a popular method to obtain new exactly solvable potentials is to apply Darboux transformation [5] 
to the ground state of an exactly solvable potential. However, when applied to the excited states, this 
transformation produces not just one isospectral potential, but a number (depending on the nodes of 
the wave function) of nearly isospectral potentials which are defined not over the whole domain, but in 
disjoint intervals Here our objective is to apply Darboux transformation to non Hermitian potentials 
and it will be shown that for such potentials, it is possible to have wave functions without nodes on the 
real line, by a reasonable choice of parameters. In the present article, we shall use this result to construct 
non trivial isospectral partners of exactly solvable complex potentials. In particular, we shall apply the 
Darboux transformation to the well known VT symmetric Scarf II potential 

V(x) = —V\sech 2 x — iVisech x tanh x, Vi > 0, V2 7^ (1) 

and generate a series of new exactly solvable non-Hermitian potentials with real spectrum. 

We note that in the case of Hermitian quantum mechanics Darboux transformation is equivalent to 
supersymmetry |10| . However this is not so in the non Hermitian case. So it is natural to ask whether 
there exists any symmetry which relates the two Hamiltonians i.e., the original and the one obtained by 
Darboux transformation. The answer to this question is in the affirmative and it will be shown that the 
two Hamiltonians are related by pseudo supersymmetry In other words the hamiltonians obtained 
by intertwining are pseudo supersymmetric partners. 

Finally we shall examine the problem of applying Darboux transformation to models with sponta- 
neously broken 'PTsymmetry. It is known |12j that models with spontaneously broken 'PTsymmetry 
exhibit a complex spectrum and all the energy levels appear as complex conjugate pairs. It will be shown 
that if Darboux transformation is applied to such a system one gets a potential with complex energy 
eigenvalues but as singlets. 

The organisation of the paper is as follows : in section 2 we briefly present the Darboux construction; in 
section 3 we construct new VT symmetric potentials; in section 4 we show that the partner Hamiltonians 
are connected by pseudo supersymmetry; in section 5 we examine the nature of the spectrum obtained 
by applying Darboux transformation to a potential with spontaneously broken VT symmetry and finally 
section 6 is devoted to a discussion. 



2 Darboux Transformation 

To make the article self contained we start with a brief review of the Darboux transformation [HUHj • A 
particle moving in the potential v(x) (real or complex) is characterised by the Hamiltonian 

* = (2) 
(The units used are h = 2m = 1 for convenience). 

If the particle is in the m th state, (i.e., m is the quantum number equal to the number of nodes of the 
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m th eigenfunction i/j m (x) of the starting potential v(x)), and the energy scale is adjusted so that the m th 
energy eigenvalue is exactly zero (E m = 0), then the Schrodinger equation reads 

H^ m = {-jL + v(x)\ - (3) 



Equation © has a potential 



v(x) = (4) 



which is regular everywhere, so that the Hamiltonian in Q may be represented as 



dx 2 ip. 

Thus if the general solution tp = ip(x) of the Schrodinger equation 

f^ + [e-v(x)]iP = (6) 

is known for all values of e, and for a particular value of e = E m , the particular solution is ip m , then the 
general solution of the equation 
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Choosing different ^> m , one obtains a series of non-trivial partners u(x) given by iJBJ, which contain all 
the states of the original potential v(x) except the m th state, i.e. the one corresponding to the eigenstate 
ip m . It is easy to observe that the partners u(x) and v(x) are related by 

v(x) = W m (x) 2 -W^(x) (11) 

u{x) = W m {xf + W' m {x) (12) 

where 

W M = (- |f) (13) 

Though ip' m 7^ at the nodes Xj, j = 1,2, 3, ■ • ■ of ip m , W m (x) has singularities at these points. However, 
since the second derivative ip!^ also vanishes at the nodes, from Q v(x) is regular everywhere. In case 
of Hermitian SUSY QM, W m is the superpotential, and W^(x) ± W^x) are called SUSY-m partner 
potentials 0. However, one can construct the partners for m = only, as for non-zero m, W m becomes 
singular and as a consequence such potentials are not defined over R but over disjoint intervals, the 
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number of intervals depending on the value of m. Thus, for m = 1, there are two potentials, each of 
them defined on a semi-infinite domain, for m = 2 there is one potential on a finite domain between 
nodes Xi and x 2 , and two potentials on the two semi-infinite domains (— 00, xi] and [x2,+oo), and so 
on 9 . In the case of non-Hermitian quantum mecahnics, W m (x) is no longer singular (and it is not the 
superpotential anymore). As a consequence the new potentials do not have singularities on the real axis 
and are defined on (— 00,00). Thus, if one of the partner potentials is exactly solvable, this formalism 
enables one to construct an infinite number of exactly solvable, non-trivial partners defined on the entire 
real line (—00, +00), unlike in the case of Hermitian quantum mechanics. 



3 New exactly solvable VT symmetric potentials 

In this section we shall construct isospectral partners of the complexified Scarf II potential. This potential 
is given by 

V(x) = -Visech 2 x - iV 2 sech x tank x, V\ > 0, V 2 ^ (14) 

and it has been studied by various authors as it is not only invariant under VT symmetry, but also 
P-pseudo Hermitian. This exactly solvable model has certain interesting properties. It has a discrete 
spectrum that admits both real as well as complex conjugate energies, depending on the relative strengths 
of its parameters V± and V 2 . The normalized wave functions for this potential are well known, being given 

by pa 

nW (2 - 2p) 

where P"' 13 are the Jacobi polynomials given by ^3] 

P^(i sink x) = _^±^±iL F (-n, n + a + (3 + 1; a + 1; z) (16) 
r(n + l)r(a + 1) 



and 



1 — i sinh 



(17) 



P = -\±\J\ + V 1+ V 2 = -\±\ (18) 



11/1 Is 

^ = -4 ± 2Vl + ^- y2 = -I ± 2 (19) 
t and s are defined with only the positive sign in the discriminant in p and q. The energy eigenvalues are 
obtained as 



E n = -(n-p-qf , 71 = 0,1,2,... < ( - ' - ) (20) 

Since V\ > 0, two cases arise for real V 2 , depending on the relative strengths of the real and imaginary 
parts of the potential : 

1. For \V 2 \ <V x + \ 

In this case the potential and the wave functions are VT invariant, p and q are real and one gets a real 
bound state spectrum. In addition to the potential (1), the wave functions ip n (x) given in (15) are also 
VT invariant. Note that due to normalisation requirements only the values with the positive sign are 
allowed in JTHJ and itl^l . 



t-1 
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2. For \V a \ > Vi + \ 

VT symmetry is spontaneously broken, as though the potential is VT invariant, the wave functions are 
no longer so. Either p or q is complex, and all energies occur as complex conjugate pairs. Real energies 
are conspicuous by their absence. 

For purely imaginary V2, however, the potential is real, possessing only real energies. 

However, when V2 has both real and imaginary parts, the potential Q loses its VT invariance, and so 
will not be considered in this study. 

We now consider the first case when VT symmetry is unbroken. Now, using the explicit solution for 
the normalized wave function (|15p. we obtain 

w m ( x) = -fM 

mi F(-m+l,i+l;c+l; Z ) [ZL> 

— (P + q)tanh x — i{p — q)sech x H ; r 

c F (— m,o;c; z) 

where b and c stand for 

b=-2p-2q + 2 (22) 

c=-2p+\ (23) 

The exactly solvable potential U^ m \x), which is isospectral to the Scarf II potential (except for the m th 
state), is obtained from the formula 

U^\x)=Wl+W' m -!5 m (24) 

In writing the last expression we have made use of the fact that if v(x) and u(x) are isospectral, so are 
V(x) and U m (x), given by 

V(x) = {v(x) - f3 m } (25) 
U m (x) ={u(x)-(3 m } (26) 
For the Scarf II potential, /3 m is calculated to be 

(3 m = (p + qf - 2m (p + q)+ m 2 (27) 

Thus this approach yields new interesting potentials, with eigenfunctions for this particular case being 
given by (see ea. <jTU|) 

In the above P n stands for P"' 13 (i sink x) and P' n denotes differentiation of P"< @(i sink x) with respect 
to x. Since P"> P is well defined on the entire real line, so also is 4> n {x). 
Let us analyse three low lying cases m = 0, 1, 2. 

For m = 

U {a \x) = -{2 (p 2 + q 2 ) - (p + q)} sech 2 x - i(p - q) [2(p + q) - 1] sech x tank x (29) 
with eigenenergies 

E n = -(n + l-p-q) 2 , n = 0,l,2,3,--- (30) 
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, 1 — i sink x\ ( p 2 ^ ( T + i sink x \ ^ q 2 ' . . 

0o = N [ (31) 



and the ground state 



Thus for m = 0, the partners belong to the family of the so-called satellite potentials. I(29|) is also a Scarf 
II potential, with a different set of parameters, and shares all the energies of QJ except for the ground 
state of V(x), which is missing in (|29|l . So, this is analogous to the Hermitian case. 

m = 1 gives the first of the non-trivial potentials. 

U {1 \x) =-{2 (p 2 + q 2 ) - (p+ q)} sech 2 x - i(p - q) [2{p + q) - 1] sech x tank x 

+ 2 (£) 2 -fc^-2 (32) 



with 



h{x) = F(-l,-p-q-X;2p+^;z 

— (p — (?) + - (— 2p — 2g + 1) sin/z x 



(33) 



The ground state of the partner 132(1 is obtained from 1(1 0(1 as 

(— p — q) + i [ ^ — p — q) sink x 

with energy 



E = - (-p - qf (35) 

It can be shown that the state corresponding to tpi is excluded from the spectrum as it turns out to be 
non-normalizable. All other states share identical energies with the original potential l|T]l. The excited 
states are obtained from (|rU|) as 



P n+2 - Pj P n+2 



V'o (36) 



with energies 

E n = - (n + 2 - p - q) 2 , n = 0, 1, 2, 3, ■ • ■ (37) 

As P"' P(i sink x) has no zeroes on the real line, so the eigenfunctions are well defined. It is easily 
observed from 1(36(1 that <f> n are normalizable. Moreover, the potential U^(x) so constructed has no 
singularity on the real line, and hence is defined on the entire domain (—00, + 00). Also, for real values 
of the parameters p and q (corresponding to real energies) the new potential, too, is invariant under VT 
transformation. 

In an analogous way, the isospectral partner for m = 2 is found to be 

U {2) {x) = - {2 (p 2 + q 2 ) - (p + q)} sech 2 x - i(p - q) [2{p + q) - 1] sech x tanh x 

+ 2 i ^\ 1 <? - 6{p - q)i sinh x Q (38) 



hj h 

with 

f 2 (x)=F(-2,-p-q-X-2p+^;z) f.'SO) 
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-*- 2 ' + 2 h (p - ,f - I fi H** I) + ■ (3 + 2„ + „) (3 - 2„ - 2,)1 (40) 



(-2p+i)(-2p + |) { " (-2p-2, + 2) 2 

For a visual representation as well as for comparison, we have plotted the real and imaginary parts of 
the potentials V , (i = 0, 1,2). In fig 1, we have plotted the real parts of V , (i = 0, 1, 2) for 

the parameter values V\ — 24, V2 = 18 while in fig 2, we have plotted the imaginary parts of the same 
potentials for the same values of V\ and V2 ■ 

Since the Scarf II potential is always VT symmetric, so are its partners [/( m ) (x), for real values of the 
parameters p and q ( V\ + 1/4 > |V^| ), as the functions f m (x) remain invariant under VT . Moreover, all 
the new potentials so constructed are defined over the entire domain (— 00, +00), admit real bound state 
spectrum, and possess all the energies of the original potential except for the m th state, if one starts with 
the m th order eigenfunction. Though m — gives the usual shape-invariant form, highly non-trivial, 
non-shape-invariant potentials are obtained for non-zero m. 



4 Pseudo supersymmetry and intertwining 

It is well known that in the case of Hcrmitian quantum mechanics, Darboux transformation and su- 
persymmetric quantum mechanics are equivalent. Although this is not so in the case of non Hermitian 
quantum mechanics, Darboux transformation can still be implemented in terms of intertwining operators. 
To see this we consider the intertwining operators A and B : 

A=l + W m (41) 

B = -^- + W m (42) 
where W m is defined by (|13fl . then the partner Hamiltonians 

H ± = -^ + Wl±W' m (43) 
can be written as £/_ = BA and H + = AB, where 

(44) 

tf + =(--*+„(*)) (45) 

Evidently, if ip n is an eigenfunction of H- with energy eigenvalue E~ , then cf> n = Aip n is also an eigen- 
function of H + with the same eigenvalue E~, escept for n — m, since in this case At{> m = 0. 

H + A^ n = (AB)Atp n = A(H^ip n ) = E~ {A^ n ) (46) 

For Hermitian Hamiltonians, A and B are mutually adjoint operators (B = A>), giving the well known 
results of supersymmetry, viz., AH- = H+A or H-A* = A*H+. 

To extend the idea of supersymmetry to include non Hermitian hamiltonians, we assume the existence 
of a linear, invertible, Hermitian operator 77, such that |11| 

B = A* = tT 1 ^ (47) 



7 



This allows one to rewrite the partner Hamiltonians as 

H+ = B*B H- = BB* (48) 

so that 

BH+ =H-B H+B* = B*H_ (49) 

From l|49fl it is clear that B maps eigenfunctions of H + to those of if_, and A(= B#) does the converse. 
Thus the mutually adjoint operators A and of conventional supersymmetric quantum mechanics are 
replaced by their pseudo supersymmetric counterparts A and B when the potential is non Hcrmitian. 
However, we would like to point out that the choice of r\ is not unique. 
To determine a form of r), let us note that a simple representation is given by 

r) = V , Vf{x) = f{-x) (50) 

It follows that for real potentials, (|50|l leads to B = A', thus reproducing the conventional result of 
supersymmetry. We note that the above results are quite general since they do not depend on a specific 
W m . We thus conclude that in the complex case the hamiltonians H± are pseudo supersymmetric partners 
of each other. Finally to cast the above results in a formal pseudo supersymmetric form let us define the 
pseudo supercharges Q and Q# in the following way : 

«-(S 0) ■ «•-*-'«*■>- (I 0) < 51 > 

Thus the pseudo supercharges Q and Q# are nilpotent and they satisfy the following closed algebra : 

H = {Q,Q#}=( H + ° ) = £ A ^j , [Q,H] = [Q*,H]=0 (52) 

We thus conclude that in the non Hermitian case the hamiltonians obtained by intertwining are pseudo 
supersymmetric partners. 



5 A model with spontaneously broken VT symmetry 

As mentioned earlier, there are two cases for the VT symmetric Scarf II potential, depending on the 
relative strengths of V\ and V2, viz., 

(i) |^2|<^i + i: 

In this case, where the spectrum is real and discrete, each state is a singlet. This has already been 
discussed in Section 3. 

(ii) \V 2 \>V 1 + \: 

In this case VT symmetry is spontaneously broken. We can choose p to be real, taking a single value 
with only the positive sign in l|18|l while q can take either of the following values 

q ± = -\±i\ (53) 

where 



V 2 - V 1 - \ (54) 



giving rise to complex conjugate pairs of energies 

E± = -nl±i[i n s , n= 0,1,2, 



1 (55) 

= Tl - p + 4 
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This case makes quite an interesting study and we shall investigate it further. Though the original 
potential (14) is still VT invariant, the partners are no longer so. We note that while the original potential 
V(x) does not depend explicitly on the parameters p and q, the partner potential is explicitly dependent 
on these parameters. As a consequence there are two partner potentials corresponding to V(x) (this is 
due to the fact that in this case both the values of q are allowed). If the Darboux transformation is carried 
out by the m th eigenstate (ipm) °f the original potential, then straightforward calculations show that 
the corresponding state will be missing in the partner. Furthermore, the entire positive (negative) sector 
E£ = — (ft* — + ifi n s (E~ = — (jjl^ — ^ — i/J>n s ) will be absent in the partner. Thus the 
spectrum of the partner potential is quite different from that of the Scarf II potential, as the former has 
only singlet complex energies. So, if one starts with the ground state eigenfunction ip^ , then the partner 
potential l|29l) is of the form 

s 2 3 s 2 3 

U°(x) = -{(2p 2 -p h -) + is} sechrx - {i(2p 2 - p H ) + s} sechx tanhx (56) 

2 8 2 8 

while the wave functions and the corresponding energy values are given by 

<t>n( x ) = — z~ p (z*)-i + ^F(a + 1,6+ l;c+ l;z) , E~ = - ifa+is , n = 0,1,2,- ■■ (57) 

c 

where 

3 1 

a = -(n + 1) , b = (n + - - 2p) + is , c = -2p + - (58) 



6 Discussion 



In this article we have explored the idea of applying Darboux transformation to non Hermitian quan- 
tum mechanical systems. In particular we have obtained a series of new exactly solvable VT symmetric 
potentials with real bound state spectra. The symmetry aspect of the potentials has also been investi- 
gated and it has been shown that they are pseudo supersymmetric. On the other hand when Darboux 
transformation is applied to a system with spontaneously broken VT symmetry the resulting potential 
admits a single tower of energy. We would like to point out that this is not a characteristic feature of 
only models with spontaneously broken VT symmetry but in general the same result would be obtained 
whenever the original model posseses doublet of states (T3) . 
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Figure Captions 

Fig. 1 : Graph of real parts of V (solid), (dotted), (small dash), (large dash). 

Fig. 2 : Graph of imaginary parts of V (solid), (dotted), (small dash), (large 
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